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SYNOPSIS 


This is the final report for NASA Grant NCC 2-711. During the final 
reporting period (June 1992 - December 1992) analyses of the ‘ . an ^ 

lateral flying qualities were made for propulsive-only flight control (POFC) o 
a Boeing 720 aircraft model. Performance resulting from compensators 
developed using Quantitative Feedback Theory (QFT) is documented and 
analyzed This report is a first draft of a thesis to be presented by graduate 
student Hwei-Lan Chou. The final thesis will be presented to NASA when it is 

completed later this year. 

The latest landing metrics related to bandwidth criteria and based on 
the Neal-Smith approach to flying qualities prediction were used in developing 
performance criteria for the controllers. The compensator designs were 
tested on the NASA simulator and exhibited adequate performance for piloted 
fliqht There was no significant impact of Quantitative Feedback Theory on 
performance of the propulsive-only flight controllers in either the longitudinal 
or lateral modes of flight. This was attributed to the physical limits of 
thrust available and the engine rate of response, both of which severely 
limited the available bandwidth of the closed-loop system. 
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Synopsis 

Flight control using throttles only to achieve safe landing for a large jet transport 
airplane, Boeing 720, was investigated. Through throttle manipuUhons en^ne 
can be used for emergency flight control for multi-engine aircraft NASA Dryden 
studiedThe use of throttles fofemergency flight control and found m genera manual 
fly-by-throttle is extremely difficult, with landing almost impossible, but control 
augmentation makes runway landings feasible. 

An auemented control developed in previous simulation study for Boeing 720(B- 
720 ) ThrottleSnly Flight control(TOFC) has obtained generally good pilot rating, bul 
whh iS t Dutch-roll damping and low control bandwidth that needs to be improved A 
controldesign technique. Quantitative Feedback Theory(QFT), and a 1, near, zed B-720 
miX wereTmployed to analytically study the B-720 TOFC and to opt.rn.ze the 
previously developed augmented control. 

The fidelity of the linear model was first examined and then modified. Poor 
open-loop characteristics were observed. The response frequency for pitch control has 
been Sased and the Dutch-roll damping doubled To -provcyhe cont^bandvndth 
substantially proved very difficult. The pitch control designed by QFT showed 
consistently^ better results for robustness and for response under no turbulence and 
turbulence 7 The roll control by QFT performed better than the roll contro built by 
stouttion st^when no turbulence was present, but not as well under intermediate 
turbulence Applying differential throttle control to engines mounted at different x-y 
planesfrom c to produce more instantaneous pitching moment may be a effective 
way to increase' 'the control bandwidth. Handling qualities of the augmented control for 
cruise and landing will be evaluated by piloted simulation fligh . 

Notation 
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Throttle-Only Flight Control 
Quantitative Feedback Theory 
nondimensional velocity-pitch coupling derivative 
nondimensional yaw-roll coupling derivative 
pitch rate (deg /sec) 
flight path angle (deg) 
angle of sideslip (deg) 
bank angle (deg) 
pitch angle (deg) 
thrust (lbs) 
natural frequency 
damping ratio 

stick input(full deflection=l unit) 
pitch rate feedback gain 


K* 

<C‘ 

(a) 

[z, w n ] 
D.R. 

s.p. 

c -g- 


flight path angle feedback gain 
sideslip angle feedback gain 
bank angle feedback gain 

transfer functons 
short form of (s+a) 
short form for S 2 + 2£ai n S 
Dutch-roll 
short peroid 
center of gravity 


+ (O 
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1. Introduction 

Through throttle manipulations, engine thrust has been found useful in 
providing some controllability for multiengine aircraft m emergency situations with ^ 
severe or complete flight control system failures. NASA has studied the use of throttl 
for emergency flight control for a range of airplanes 1 ' 6 - Many multiengine airp anes 
exhibited some degree of useful control capability with the throttles. In genera , 
manually flying an aircraft using throttle-only requires tremendous pilot workload and 
makes landing extremely difficult to almost impossible. Control augmentation using 
feedbacks antf direct coupling of the throttle command to stmk/thumbwhed motton 
has greatly improved the flying qualities and successful landings could be ma . 

An augmented control developed in previous studies for B-720 TOFCH 
implemented^ a high fidelity B-720 flight simulator(Fig.l &2), had obtained genera y 
good pilot rating. However, its light Dutch-roll damping and low 

desired to be improved for better handling qualities. The pnmary aim of jfois study) is to 
use an alternative control design technique. Quantitative Feedback Theory(Q T), 
linearized B-720 model to analytically study the B-720 TOFC and to improve e 
roll damping and control bandwidth for better handling qualities. 

The OFT technique 18 - 23 was chosen because of its unique features and of the 
good insights it provides throughout its design process. It allows designers to specify a 
desired close loop frequency response with a desired control bandwidth and damping 
characteristic, and then designs controllers and prefilters to meet the prescribed 
specifications. The desired performance can't always be achieved ^ thin ^ 6 ^ 
control actuation and rate limits; nevertheless, the technique provide 
show clearly what has prevented the implementation. The designer can thus hav e bet 
understanding and make better judgment for the system being investigated. 

In this report, the concepts of throttles only flight control ^ infroduced ™d an^^ 
augmented throttle-only flight path control for approach and! Ending for B-720 is bi 
using the QFT technique applied to the linear B-720 model. The fidelity of the 
linear model is first examined. The B-720 engine model and the open-loop bare airframe 
response are then analyzed and discussed. Augmented control design using QFT 



presented in a step by step procedure and Its results are compared with the one from 
the previous simulation study. Conclusions and recommendations are then made for 
future studies. 


2. Concepts of Throttle-Only Flight Control 

The propulsion system of a multi-engine aircraft can be used for heading and 
flight path control. Differential throttle is applied to control roll through yaw and 
symmetric throttle to control pitch. Speed control by throttles becomes ineffective when 
control system fails. Other means may be used to change airplane's speed. Throttles are 
coupled to stick/ thumb wheel for easier and more conventional control handling. 

Pitch Control Symmetric throttle induces a speed change which in turn generates a 
pitching moment change through speed stability effects-C mu , this is the primary source 
of pitch control. A smaller but immediate pitching moment change is contributed by the 
offset of thrust vector to the c.g. of the aircraft. Engines mounted at different x-y planes 
from c.g. can produce additional instantaneous pitching moment by applying 
differential throttle control. This was proposed initially as a control strategy, but has not 

been investigated yet. 

Yaw-Roll Control Differential thrust generates sideslip, which in turn generates 
rolling moment changes through wing dihedral and sweep effect-C //3 . Roll is controlled 
by applying differential throttles to obtain desired bank angle to make turns and 
heading changes. 

Speed Control Retrimming speed by throttles becomes ineffective with primary 
control surfaces locked due to control systems failure. When control system failure 
occurs at speed other than landing speed, retrimming to an acceptable landing speed 
may use other techniques, such as moving c.g. aft, lowering flaps, extending landing 
gears and variable stabilizer control to reduce speed. 

Couple Throttle Command to Stick/Thumbwheel Motion Direct coupling of the 
throttle command to stick/thumbwheel motion has eased the pilot's handling of control 
and allowed him to control the airplane in a conventional fashion, i.e. pitch up with 
stick forward, pitch down with stick aft, and etc. 

3. B-720 Linear Model 

The B-720 linear model was derived from perturbations of the full non-linear 
equations of motion about trim and was completely decoupled in longitudinal and 
lateral dynamics. 

The fidelity of the linear model was examined by comparing the open-loop 
response of the linear model with the nonlinear as shown in Fig. 3 & 4. In Fig. 3, the 
longitudinal response of the linear model was about 30% less in magnitude than the 



nonlinear's, and was, therefore, modified by a correction factor of 1.3. The result after 
modification is shown in Fig. 5. The linear design used a computer control package 
"CC" to assist the design. 

The lateral response in Fig. 4 showed the linear model would closely follow the 
nonlinear model as long as the small perturbation assumption was not violated, i.e. the 
command input should be of small magnitude and short duration. 

Fig. 3 & 4 also showed that for the nonlinear model, a flight path angle command 
would induce little coupling in roll /yaw while a bank angle command induced 
pronounced pitch coupling. Coupling between longitudinal and lateral modes was 
completely absent for the linear model. 

4. Engines and Bare Airframe System Analysis 


4.1 Engines 

The spool-up and spool-down engine dynamics for the B-720 engine are shown 
in Fig. 6a. The empirical transfer function developed is given in short form notation by 

QZ(lb6) = 275 

(,55)(5) 

The above equation is illustrated in Fig. 6b over low frequency ranges up to 1 .0 
rad /sec. 

4.2 Bare Airframe 

It is apparent from the engine bode diagram that severe bandwidth attenuation 
occurs beyond frequencies of 1 rad /sec. It may not be possible, therefore, to increase the 
closed-loop bandwidth beyond 1 rad /sec within the range of available thrust. 

This can be seen in the pitch rate "q" to thrust "z" transfer function of the bare 
airframe shown in Fig. 7. The full-order transfer function G^ /KC) shows that 80 db of 

gain must be added to yield a crossover frequency beyond 1 rad /sec. This corresponds 
to 10,000 lbs of full thrust from each engine, which would not be practical for approach 
and landing. 

A low order fit to G? de8/sec) is also depicted in Fig. 7 and is very accurate near the 

phugoid frequency. Piloted flight of the unaugmented aircraft was consistently level 3. 
The main difficulties were the lightly damped phugoid and the low bandwidth throttle 
control. 

The accuracy of the low order fit near the phugoid frequency means that, to a 
first order approximation, the phugoid frequency and damping are found from. 



2£o> n = -X u + 


M„(X a -g) 


w„ 2 = 


-g(Z u -^ZJ 
Mq 

U„ 


and for conventional transport aircraft can be shown to be roughly proportional to M u . 


It should be strongly noted here for the classic case of M u =0 and for negative 
values of M u (Mach tuck) that the aircraft cannot be practically flown with throttle alone 
unless rotational control in pitch is added. Difficulties will also be encountered as M a 
becomes small (aft c.g. location.) Both of these cases require the addition of an effective 
rotational controller about the pitch axis. This may be achieved using differential 
inboard and outboard thrust, provided the inboard engines are a different distance from 
the aircraft xy-plane than the outboard engines. These configuration charactenstics 
determine the innate capability for throttles-only piloted control. 


5. QFT Control Design 


An overview of QFT design technique is presented in Appendix B. 

In the following, a brief procedure of QFT technique is introduced. QFT is a 
frequency domain control technique that uses a fairly direct and transparent design 
approach. In order to apply QFT, a system must be modeled as a SISO/MISO unit 
feedback system. For M1MO systems, a m x m MIMO system can be converted into m - 
equivalent SISO/MISO system. QFT technique allows designers to specify a desired 
performance specification and then incorporates the plant uncertainty and system 
disturbance to form design constraints - performance bounds and a U contour. 

After the design constraints are formed, they are displayed on the Nichols Chart 
together with the plant transfer function. A controller will then be selected to reshape 
the plant transfer function to satisfy these constraints. After satisfying all the design 
constraints, if possible, the system is guaranteed robust over the full range of plant 
uncertainty but may not completely meet the performance specification yet. A prefilter 
is usually required to reposition the compensated system, shown on a frequency plot, to 
fully meet the specification. 

A QFT computer control package was used to assist the design, the program can 
apply QFT on minimum-phase plants only, i.e. the plants should have no zeros in the 
right half s-plane, and only the gain of the desired close-loop performance will be 
specified and then satisfied. For nonminimun-phase plant, the phase of the desired 
close-loop performance shall also be specified and satisfied. 

The QFT technique for B-720 TOFC design is presented in steps as follows: 



5.1 Syste m Modeling 

TO apply QFT, a system should he in 
The B-720 pilot simulation block d.agram of flight path ana name g 

presented below: 






form as shown below: 
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Bank An f 1e Control Block Di agram wit h Inner Loop Clo sed 



52 Performance Spec ification 

OFT allows designers to specify a desired close loop frequency 

responsewith an upper^ound, B u . a lower bound, B. and a tolerance, d and 1 a 

maximum peak magnitude, M m (see Fig. B4). The tolerance, 6 ts s P ec * e ^'° obtam 
robust performance, and the maximum M m to obtain a desired system damping. 

The upper bound, B u , is usually described as an underdampec ' « s P onse of a 
second ordersystem, while the lower bound, B u an overdamped response. B u can 
determined by satisfying some or all of the desired performance i spe 
UI K ( and B L by satisfying t r , For example, if the desired 

tp,VrNmfrdu l y j o a nd tplO sec, then B u for the 

specifications are: damping ratio z=0.6, rise time tp-asec, s- 

flight path angle and bank angle feedback loops can be modeled as follows. 









n 


« n 

■ g = 0.6 , and t p = 7 = 2 

®n V 1 — 


= = 5 see 




•a,.-. 78 rad /see, M„ = 1 - 1.09 and l s . -i- = 87 sec S 10 sec 

since and 


.61 


therefore, T B> = JT + 94S+ 61 
a zero at S = -5 is added to T„ u to widen the <5 R at high frequencies 

thus B - T = 3-05 ( s + 5) an( j L m m =.30 db 

thus Bu (s 2 +.94S+.6l) 

And Bl can be modeled: 

— = 8.7 sec and g*l, 


vt s = 

£<°n 

co <s 0. 76 rad / sec 


K 


.385 


T XV ■ 

therefore, T* = — = (S+.55) (S+.7) 

where cr n should be i co n , 

and t s of T„ need to be s 10 sec 

a pole at S = -2 is added to T Bl to widen the d R at high frequencies 

0.77 

thus B t - T a - f s+ 55 ) (s+.7|fs + 2) 

The magnitude of B„ and Bu and «„ at each frequency can thus be obtained and are 
shown in Table 1. 


IdUlc 1. 1 iciiiiim 

Freq uency (rad / sec) 

zlZ r , 
0.1 

0.3 

r 

0.5 

0.7 

1.0 

2.0 

5.0 

Bu(dB) 

0.0 

0.0 

0.3 

-1.0 

-4.0 

-15.0 

-30.0 

BL(dB) 

-0.2 

-2.0 

-4.7 

-7.1 

-12.4 

-24.1 

-45.0 

<5„ (dB) 

0.2 

2 

5 “I 

6.1 

8.4 

9.1 

15.0 


The performance specification obtained above is the desired closed-loop 
response for y - and q> -loops with no pilot interface. Since these two feedback loops will 
become the pilot open-loops, additional specifications: a desired control bandwid ( 
BW) of 2 rad /sec for landing tasks for transport aircraft(see ref. 17) and a k/s gam slope 



near crossover frequency, need to be added in order to obtain good p, lot handling 
dualities The OFT technique describes here applies only to minimum phase system . 
Smtinaminimum phase system the magnitude of the frequency response completely 
specifies the transfer function, hence the BW galn , the bandwidth demined by the 

r of ^ hk ran snlplv determine the system bandwidth, the D u ana dl acquireu 

[Xreforel needed to be r^d S meet the additional specifications. The performance 
specification is finalized as shown in Table 2 and Fig. 8. 


laDiez.. rmcii pci 
Frequency(rad/sec) 

1U1 11 IC11 
0.1 

0.3 

0.5 

0.7 

1.0 

2.0 

5.0 

Bu(dB) 

17.0 

17.0 

17.3 

16.0 J 

13.0l 

2.0 

-13.0 

BL(dB) 

16.8 

15.0 

12.3 

9.9 

4.6 

-7.1 

-23.0 

i r n 

6 a (dB) 

0.2 

2 

5 

6.1 

8.4 

9.1 

15.0 


S-3 Design constraints 

The performance bounds constraint is a curve on Nichols Chart that constructed 
by mapping^ the template in s-plane, formed by plant uncertainties, onto the ic 
Char^ tif match the performance tolerance, d„(see Fig. B-2 & 3). 

constraint guarantees the variation of the system response due to plant uncertainties 
will be no greater than f>„ When the template of plant uncertainty in s -plane is 

mapped onto Nichols chart(joi-plane), the template transforms into d |^ , *^«P e f ° r 
each different frequency. Thus, there is a performance bound for each frequency. 

The U contour is a M<ircle on the Nichols Chart that has the magnitude of M„ 
With nart of the circle stretched for uncertainty at high frequencies. By having the 
r y tm not peid Seating the U contour, the system's damping will be guaranteed no less 

than the damping correlating to M m . 


5.4 Airplane Parameter Unc ertainty 

Four configurations were given to study the approach and landing of B- 
720 throttle-only flight control. The flight condition of these configurations are 
summarized in Table 3. Configuration 1 was used as the nominal confiuration 

for control design. 




Config. 

Number 

Weight 

(lbs) 

Altitude 
(Ft MSL) 

Airspeed 

(Knots) 

Flaps 

(%) 

Gear 

up /down 

j 

1 

140,000 

4,000 

' 160 

0 

U P 

2 

140,000 

4,000 

145 

30 

U E 

3 

160,000 

4,000 

175 

0 

U P 

4 

140,000 

i 4,000 

155 

30 

U P 

Table 3. Fligl 

ht condition of the four configurations given for B-720 1UFI approacn ana. 


study 


A robust controller is highly desirable for systems undergoing constant changes 
such as aircraft. QFT incorporates the plant parameter variation(uncertainty) into 
design by converting them into design constraints and then builds controllers to satisfy 
the constraints. An example of a plant with parameter variations is illustrated below. 


Example 

For a plant transfer function 

CfS'v = Ka _ where the parameter variations are: 1 < k < 10 and 1 < a < 10 

(s + a)’ 

then, 

1 „ 10 ° 

G(S) min.-— and GfS) max. - s + )()0 

For inner-loop transfer functions GJ, and G*. , the parameter variation given by the four 
configurations can be expressed in the same manner . 

For longitudinal flight path angle control: 


The G£ (d ?! ) , of the nominal configurationfconfig. 1) is: 

^In (° c 8) 


,01(. 203137, 3.01] 

G £K*> conflg ' 1 = (.562) [.624, .111] [.441, 1.57] (5.25) 


and the min. and max. Gj^f^ are: 

. .0053 (.162) [.35 , 3.0l] 

G e,„(^g) min,= ( 40) [.42 , 1.48] [.66 , .0l] (5.19) 

max .01 (-28) [-46, 3.43] 

G e in «* g ) max - (_ 58 ) j- 45 5 L57 ] [.92 , .14] (5.24) 


For lateral bank angle control: 


The G* (dcg) of the nominal configuration(config. 1) is: 

Pint**) 

. 09 [. 47 , 3.65] 

G;<X, nomunal = ( 98) [ 81 , , 1.07] (5.65) 

and the min. and max. of Gj^ g) are: 

, .06 [.45, 3.651 

G A.3» mm ' = (.98) [.60 , .15] [-24 , .93] (5.01) 

.09 [.61 , 4.33] 

Gjjdeg) max -= ( 10 3) [1.0 , .20] [.29 , 1.09] (5.02) 


The minimum and maximum values of transfer functions, and GJ^, are given 
to form the uncertainty template which will then be used to determine the performance 
bound constraints as explained in section 5.3. The QFT control package allows the 
designer to input plant parameter variations by entenng the transfer function s 
maximum and minimum - gain, first order poles and zeros, and second order pole 
zeros. However, due to the software's limitation of handling the quantity of 
uncertainties, some uncertainties have to be averaged instead of a maximum and a 
minimum. There are tradeoffs between plant parameter variations and performances. 
The wider spread of the parameter variation, the more restricted the constrain s, 
consequently more compensation is required. Therefore, the performance specifica ion 
may need to be relaxed when there is not enough control power to provide all the 
compensation that's required. 


S-S Controller design 

Poles /zeros /gain compensation may be required to reshape the inner loop transfer 
functions, G' and G* (see control block diagram in sect. 5.1), to satisfy the performance 
bounds and the U contour constraints. An added gain will raise the transfer fimction curve, 
a zero will bend the curve to the right, and a pole will bend the curve to the left, 
poles/zeros/gain that selected to achieve the reshaping forms the controller, G c . G^ and 

G* after reshaping, become, respectively, V er and L' e „ these are the open-loop transfer 
functions of the flight path angle and bank angle feedback loops (note that 
■ r * nr an H p = G^ 1 " * G* ref. to block diagrams in section 5.1.) L er and L e ^, 

should be kept on and above the Bo(jw0, for each frequency, wj on V er and L e , . ms is to 
assure robust performance. U, and L*, should also not to penetrate the U contour in order 
to' obtain a desired damping characteristic. 


For Longitudinal Flight Path Angle Control: Transfer function G 

i • i f" 


and its 


performance bounds, Bo(jw0, and U contour are displayed on a Nichols Chart in Fig. 9. 


Notice in Fig. 9, all the frequencies on G^ are below their corresponding Bo(j Wi ). 
Therefore, reshaping is required. Since pure gain compensation is always preferred if 
possible, hence, G' was first raised by increasing the gain of G c (= G?" here) to 16, until 
it touched the U contour as shown in Fig. 10. Poles/zeros compensation was then 
added to further reshape the G^ to avoid the U contour while satisfying all the Bo( ]Wi ). 
However, no poles /zeros compensation was found that would do so. A zero at .1 
rad /sec was first added, which pulled the whole G^ curve to the right of the U contour, 


however, once a pole is added, no pole at any location would prevent the G^ curve 
from penetrating the U contour, which is least desired. Since compensator with ^ordy one 
zero is phisically unimplementable. Therefore, the compensator, G c , for the flight path 
angle feedback loop, G?" , is ,a pure gain of 16 and only the U contour constraint was 

satisfied in order to obtain a desired system damping; the performance bounds, hence 
the system robustness was left unsatisfied. The frequency response of the close-loop 

transfer function,^, where ^ =L er /(I + 14 ,)= (ft *6^/(1^; *G^), is shown 
in Fig. 11. It can be seen in Fig. 11 that <5 T , the spread between T max and T min , had 
exceeded the <5 n over the frequency range .1 to .7 rad /sec as a result of leaving L er not 
satisfying the performance bounds over that frequency range. To have any frequency, 
a) on L r e higher than its corresponding Bo(jwi) will result in d T (j<u) ^ while 

lower than Bo(jw s ) result in d,(jo>) £ tf„(ja)). It can also be seen in Fig. 10 that further 
modification is required to fully meet the prescribed specification. The frequency plot of 
the close-loop response after adding a prefilter is shown in Fig. 1 . 


For Lateral Bank Angle Control: Transfer function G* n and its performance 

bounds, Bo(jwi), and U contour are displayed on a Nichols Chart in Fig. 13. Notice that 
GJ is not only below all the Bo(jwi) but also penetrate the U contour. Therefore, more 

than just a pure gain is required to reshape G^ . A controller, Gft =(s+.15)/(s+1.5), was 
added to G* to reshape it and had avoided it from penetrating the U contour but not 
successful in satisfying all the Bo(jwi). After reshaping, L' e , is shown on a Nichols Chart 
in Fig. 14 and the frequency plot of the close-loop transfer function, T #in , where 
T* _ |* e / (i + 1* ) = (Gft * GJ ) / <6ft * Gj n ), is shown in Fig. 15. Again, prefilter is 
required to further meet the prescribed specification. 


5-6 Prefilter Design 

Gain /poles /zeros compensation may also be used for prefilter to reshape 
the close-loop frequency response, T r r |n and Tj n . It can be seen in Fig. 11, T rin 


needs both gain and lead compensation to raise the whole curve and to 

gain of 6 proved most effective in increasing the bandwidth and to meet t 
prescribed specification. 1J, after adding the prefilter is shown in Fig. 12. 

Fig. 15 shows the close-loop frequency response of T£ after reshaping but 

with no prefilter applied yet. Sufficient gain is available here. A lead 

compensator of (S+l) /(S+2) is added to haunch up the severely detenora ed 

curve at frequency over 1 rad /sec, and a lag compensator of (S+.25) / (S+.15) is 

added to steeper the gain curve at low frequencies as to prov.de a smoother k s 
addedto steep g ^ dose-loop response after adding the 


curve iui u i , , j • 

prefilter is shown in Fig. 16 and the prefilter selected 1 

15(S+.25)(S+l)/((S+.15)(S+2)). 

6. Results and Discussion 

The objective of this study is to improve the B-720 au ^ ent ^ h “^ h . roll 
rnntrol thus the primary task is to increase the low bandwidth and the light Dutch ro 

and/or different control theories may change the response. 

For bank angle control, (1 feedback was found most effective in increasing _Dutch- 

controls are presented in section 5.1. Table 4 lists the companson of the natural 
frequencyfwj and damping ratio© of the dynamic modes of the bare airframe, an 
theaugmented control built by the simulation study and by QFT. 


longitudinal m odes 
Bare a/f 

Simulation augmented control 
QFT augmented control 


density 

(1.4e-6) 

(4.7e-6) 

(3.4e-6) 


phugoid 
1.04, .13] 
[.52, .24} 
[.62, .32] 


s.p. 

[.65,1.4] 
[.52, 1.5} 
[.46,1.6] 


engines 
(.55) (5.2) 
(.4) (5.2) 
(.3) (5.2) 


Lateral modes 
Bare a/p 

Simulation augmented control 
QFT augmented control 


spiral 
(1.1 e^t) 
[. 73 , .35] 
(.39) 


engine 

(.55) 

(.45) 


D.R. 

roll 

engine 

[.12,. 99] 

(1) 

(5) 

[.15, .99] 

(1) 

(5) 

[.29,1.0] 

(1.5) 

(5) 


Table 4 Natural frequency and damping ratio of the dynamic modes of the b ag 

airframe and of the augmented control by s imulation atud y and by QF i 

For longitudinal control, pure gain compensation was used. Since the short 
period mode has frequency around 1.5 rad /sec, which is beyond the frequency that 
throttles can control, therefore, short period damping was traded for phugoid damping 
and frequency. The short period damping decreased from .52 to .46 rad /sec while the 
phugoid damping and frequency increased from .52 to .62 and from .24 rad / sec to 32 
rad /sec, respectively. This increase of response frequency can also be depicted from 
flight path angle response in Fig. 17. 

For lateral control, pole /zero compensation was used. The Dutch-roll damping 
was almost doubled, from .15 to .29. The previously build lateral phugoid mode, [.73, 
.35], which combines the spiral and the slow engine mode, was replaced by two real 
root modes, (.39) and (.45), both with higher frequency. 

All plots from Fig. 17 to 22 were obtained from simulation runs at approach and 
landing conditions with major control surfaces(ailerons, elevator sand rudders) locke 
except the electrical, the mechanical systems and the landing gear remained operative. 

The response of the augmented controls built by the simulation study and by 
QFT were compared in Fig. 17 and Fig. 18. The Dutch-roll oscillation was removed after 
the damping was increased by QFT(see Fig. 18), while the response frequency increased 
a little, from 0.35 rad /sec to .39 rad /sec. In both figures, the control magnitude and rate 
has been increased but remained reasonable. 

Turbulence Response The response of the flight path control under 
intermediate turbulence is presented in Fig. 19. Because of the gust randomness, more 
than one simulation run was made to examine the tracking integrity under turbulence. 
The presently build control showed consistently better results. 

In Fig. 20, the bank angle tracking of the QFT augmented control performed not 
very well under intermediate turbulence. Each gust caused too much tracking 
deviation. It was thought that the not well-behaved tracking is because of the larger Kp 
gain being used, the Kp was increased from 1 to 4 as can be seen from the list of the 
compensation built by the simulation study and by QFT: 



For longitudinal flight path angle augmented control: 


Simulation studv compensation: 

PF-10 , G =1 , 

T 

G® 1 ' = 10 

K r =1 , 

K q = 4 

QFT compensation: 


PF = 6 , G®" = 16 , 

Ge T ;=i 

K =1, 

o 

kO 

ii 

O’ 


For lateral bank angle augmented control: 

cumulation Study compensation : 

PF-40, B£-l, G' 1 ; - 1 , G£ =1 

K, =1. 

OFT compensation : 

PF = 2.5(S+. 25) / (S + 1 . 25) , G*" = 1 


K p =-5 


K, = (S+.15)/ (S + 1.5) , 


Ge U = 1 , 

K„ =4, 


This larger Kb multiplies the sideslip from gust plus the sideslip output from the 
airplane^ four lime larger before it was feedback to the airplane. This has dramatic effect 
onto bank angle output due to the fact that 0.2 degree of sideslip angle would Unduced 
approximately 10 degrees of bank angle owing to the large Q b of B '^°- S ^ nc ® fo 
TOFC b feedback is the only parameter that can effectively increase Dutc - 
damping, a compromise seems necessary between lateral phugoid tracking and Dutch- 

roll damping. 

However, during investigation of resolving the above mentioned bank angle 
tracking problem, it was found that the p being feedback in the simulation was the P at 
the c.g. instead of the p at the nose boom. The nose boom p is actually measured an 
feedback in a real airplane. The nose boom p was then modeled into the B-720 
simulator, which improved the bank angle tracking under turbulence. The p at nose 

boom has two extra terms than the c.g. p, one-a function of: roll rate 

of yaw rate! It was thought that the extra yaw rate term might have s abi hzed the bank 

angle tracking (the roll rate term was too small and thus being neglect .) 

seems contradictory to the preliminary investigation, which found yaw rate feedback 

does not help either Dutch-roll or lateral phugoid mode. Adding yaw rate feedback to 

bank angle control will be investigated after this preliminary report. 

Fig. 21 & 22 show the system response to configuration variations, for flight path 
control and for bank angle control, respectively. The robustness of the flight path 
control has been improved by QFT as shown in Fig. 21. Fig. 22 shows the Dutch-ro 

oscillation in simulation compensation was taken out by QFT com P®^ 1 ^ d ^^ r r/ 
the tracking became worse. After this preliminary report, using yaw rate feedback 



using Cher types of compensation will be investigated to attempt to improve the bank 
angle tracking problem. 

7. Conclusions and Recommendations 

impossible, but with control augmentation runway landing is feasible. 

• , n 1oe^72 C 0 On wa°' mVesUgate ^entedlto^ 

SffiK ^“muUHon stud 8 y has maTe successful s-ulat.o^anding, however, 

QmechntquT has improved the 

bandwidth and the Dutch-roll damping. 

For pitch control, “££5^^ *“ 

phugoid tracking under intermediate mode, 

adequate compromise is required ew , trarkinp- for bank angle control under 

Further investigation will be made on ba ^ang i mD rovethe bandwidth substantially 
turbulence. The bare airframe analysis s o P frequency stability effect 

would be difficult due to the engine performances and y 

that provides pitch and roll control. A strategy ^P^^ ounte( j en gi n es to produce more 
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Figure 1. Boeing-720 
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Figure 4. Comparison of the lateral open-loop response of the B-720 linear 
and nonlinear model, 5% differential throttle, nominal configuration 





a) Engine spooling block diagrams b) Engine 6 to 6j c 


Figure 6. Engine spooling block diagrams and Engine 6 to <5 Tc 
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Figure 7. Pitch rate to thrust bode 



, .. . , . b) Performance specification in frequency domain 

a) Performance specification in hme domain ' r ^ 


Figure 8. Preformance specification in time- and in frequency-domain, for flight path 
angle and bank angle controls. 
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Figure 9. Transfer function , its performance bounds B(jw), and U contour o 
Nichols Chart 









Figure 10. Open-loop transfer function, \J 0 ^, on Nichols Chart 



























Figure 13. Transfer function Gj |#/ its performance bounds B(jw), and U contour on 
Nichols Chart 






Figure 14. Open-loop transfer function, L^, on Nichols Chart 
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Figure 16. Frequency plot of the close-loop transfer function T* n with prefilter 
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Figure 17. B-720 augmented flight path angle control, step flight path angle response 
with no turbulence, nominal confiuration 
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Figure 18. B-720 augmented bank angle control, step bank angle response with no 
turbulence, nominal confiuration 
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Figure 19. B-720 augmented flight path angle control, step flight path angle response 
with turbulence, nominal confiuration 
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Figure 20. B-720 augmented bank angle control, step bank angle response with 
turbulence, nominal confiuration 
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a) Simulation compenstion 



b) QFT compensation 


Figure 21. Robustness of the fligh path angle control 





a) Simulation compenstion 



Figure 22. Robustness of the bank angle control 
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The B-720 piloted simulation can be represented by the following block diagram: 


Flightpaih Angle Control Simulation 



where PF=1, K" =16, K ' =1 ' =1 ' Kq=60 



where PF=2.5*(s+.25) / (s+1.25) ^ 

Gj"=i,K,=(s+.15)/(s+1.5), G e T ; =1, K />= 4 

The "AIRCRAFT" in the box above represents both the engine ‘. and 

The engine is approximated by a transfer function, G eng , and the bare airr ram y 

represented mathematically by a single quadruple, Pa/c, shown as follows. 

























I nngitudinal Dynamics 



-I 

AjJB] 

throttle inputs equal 

[A 

| B t (column)] 

. 

C j Dj 

= 

C 

j D = 0 j 


x = [q( deg/ sec) | a(deg) \ v(kts) | 6>(deg) { h(ft)] 

U - [Z^ , ef ( ll>S) 1 Z iDbd kf ,(lbs) | z i«bd Hgh( lbS ) 1 Sghflb 5 )} 

Uj = z(lbs) [used when all throttles have same command ] 


l .ateral Dynamics 


rxi 


[A 1 B] 


x‘ 


A B 2 

■ X ■ 


— 

-_ T -- 
C | Dj 


111 


cTd-o 

U 2 


where 


1 a/c 


A 1 B 

c[d 


...for four engine inputs, u 
.for one total engine input, u 2 


A | J^2_ 1 
= C|D = o] 

x = [p(deg/ sec) r(deg/sec) | f5( deg) <K de g)] 
y = j p(deg/ sec) ! r(deg / sec) j ^(deg) j <K de §) 

U =[ Z oUbdleft( lbS ) ! Z mbdleftO bS ) ! Z .nbd right ( lt)S ) ' Z ootbd right ( lt>S )] 

u 2 = [z(lbs)] , where z = z outbd left + z inbd i cf t + ( _z inbd right) + ( z 


'outbd right 


) 


rhe B metrix has four columns, each column is to be multiplied by the ^ 

s: i-:i— =; r.;s." *4* ,«• » ssr 

right engines are given same amount of command but in opposite direction), the B matr 
laferald^rnamics becomes another single column matrix, B.Each row element mB„s the sum 





the positive value of column 1 & 2 and negative value of column 3 & 4 of each row inB - The 
apen-loop configuration then becomes P=Pa/c*Peng, where Peng is the quadruple form of the engine 

transfer function, The quadruples for four different configurations were obtained as 

described in Reference. 


Phe flight conditions for each configurations are summarized in the table below. 

Configuration Summary 


Config. 

Number 

Weight 

(lbs) 

Altitude 
(Ft MSL) 

Airspeed 

(Knots) 

Flaps 

(%) 

Gear 

up /down 

i 

140,000 

4,000 

160 1 

0 

U P 

2 

140,000 ” 

4,000 

145 

30 

r up 

3 

160,000 

4,000 ~ 

175 1 

0 

U P 

4 

140,000 

4,000 

155 

30 

up 


The transfer functions were obtained from the quadruples using System Technology s CC 
Program. These aircraft transfer functions are listed here with each respective row of numbers 
designating the corresponding configuration transfer function values. The nominal 
configuration, number 1, is represented by values in each row 1 below. 


l ongitudinal Transfer Functions 


N^:f ec) = N^ g r ) /A 


'z(lbs) 


long 


wr(* g) - wr(deg) / a 
^z(ibs) ^z(ibs) ' ^long 

■£ 36E" 04 " (6) " C-lTl 7E - 05) ” ” (0.40)_ _ (0.61) 
Z33E-04 (6) (1.4E-06) (0.635, 0.563) 

^ z(lbs) 1. 976E-04 (0) (0.292) _ (_0.644)_ 

1.955E - 04 (0) (2.68E ' 06)~ ]o.819^ ^ ^ 0. 508) 

~i: 796 E - 05 _ ~ ~ "(0) " ~ ~ "(07 203) _ _ (0.370_ l 3. 008) _ 

de -1.819E" 05~ (0) (0.364) _(2-_255)^ 

N z(i^ = 2:"l 30E - 05~ "(0. 1 67)" " (0. 35 1 ~ " _ 3.038) 

1/47013-05 (0) (0.261) (0.460 , 3. 426) 


(L438E-05) 

(3.918E-02 , 

0.130) 

(0.652 , 

1.382) 

(L101E-05) 

(7.423E-02 , 

0.147) 

(0.596 , 

1.375) 

(3. 949E - 02 , 

0.118) 

(0.649 , 

1.301) 


(1.878E-05) 

(7.190E-02 , 

0.138) 

(0.588 , 

1.279) 


lateral Transfer Functions 


nks - nks / a, 


N*Sg> - Ngg 1 / k 


l^j^(dcg) _ 


-1.58E-03 

(-.0805) 

(.927) 

-1.59E-03 

(-.0922) 

(.904) 

-1.43E-03 

(-.0723) 

(.981) 

-1.44E-03 

(-.0879) 

(.940) 


3.19E-04 

(.468 , 3.65) 

2.15E-04 

(.611 , 4.17) 

= 2.89E - 04 

(.447 , 3.96) 

2.04E-04 

(.593 , 4.33) 


(.0001) 

(1.01) 

(.116, 

1.05) _ 

(.006) 

(1.05) 

(.067, 

0.93) 

(.0028) 

(1.06) 

_C.114, 

_L08)_ 

(.0065) 

(1.09) 

(.060 , 

, .944) 


The engine transfer function for all configurations is given in short form notation by 


275 

(a 55X5) 



Configuration Storage Table for Quadruples 

*Pxxxx.4U : Quadruple with four engine inputs 



1 

Config. 

Number 


Quadruple 

Pa/c 

Longitudinal 


1 


P1000.4U 



1 


PI 000.1 u 

Longitudinal 


2 


P2000.4U 



2 


P2000.1 U 

Longitudinal 


3 


P3000.4U 



3 


P3000.1U 

Longitudinal 


4 


P4000.4U 



4 


P4000.1U 

Lateral 


1 


P5000.4U 



1 


P5000.1 U 

Lateral 


2 


P6000.4U 



2 


P6000.1U 

Lateral 


3 


P7000.4U 



3 


P7000.1U 

Lateral 


4 


P8000.4U 



4 


P8000.1 U 



PI 00 
P200 
P300 

P400 

P500 

P600 

P700 

P800 


































Appendix B: Overview of OFT Design Method 

Quantitative Feedback Theory(QFT) is a theory that emphasizes the use of 
feedbacks to achieve a desired system performance while incorporating plant 
uncertainty and plant disturbance into control design 17 * 23 . The QFT technique provides 
several advantages: 

• QFT technique provides good insight throughout the process of controller design. 

• QFT technique allows the designer to specify a desired close-loop performance. 

• QFT technique allows the designer to incorporate plant uncertainty and plant 
disturbance into controller design. 

• QFT technique will design a robust controller that meets the prescribed performance 
specification over the full range of plant uncertainty. 

The basic design procedures of the QFT technique for minimum phase systems 

are: 

1) Model the system into a unit feedback SISO/MISO system to apply QFT. A m x m 
MIMO system can be converted into m-equi valent SISO/MISO systems and the 
coupling between loops can be considered as disturbance input. Fig. B-l shows a 
general SISO/MISO system for QFT application. 

2) Specify a desired close-loop frequency response. Fig. B-2 shows the construction of a 
desired close-loop performance, with an upper bounds, B U/ a lower bound, Bl, and a 

tolerance, 6r, and also a maximum allowable close-loop magnitude, M m . The upper 
bound is generally described by an underdamped closed loop response and the lower 
bound an overdamped closed loop response. Fig. B-3 shows the construction of a 
desired disturbance performance specification which only has a upper bound. The 
objective is to design a controller such that the variation of the response due to plant 
uncertainty lies within the specified boundaries and the effect of disturbance is 
minimized, that is to have: 

Bl(o))s |TR(joojsBu(co) 

|TDi(jco|s|MD(jcoJ 
|Td 2 (jco j s |m d (joo ^ 

3) Convert the performance tolerance, 6r, and the maximum M m onto Nichols Chart to 

form design constraints: performance bounds, B(joji), and the U contour. 

i) Performance bounds are constraints to ensure that at each frequency, the 
variation of the system response, 6t, due to plant uncertainty does not exceed 
the prescribed performance tolerance, 6r. 



ii) The U contour is a constraint to ensure the system damping will be kept no 
less than the damping correlates to the maximum M m . Figure 4 shows the 
construction of a U contour. 

4) Reshape the plant transfer function. With the plant transfer function, Po, and the 

design constraints: performance bounds, B(jcoj), and the U contour, all displayed on 
Nichols Chart, reshaping may be done by adding poles /zeros /gain compensation to 
Po to reshape Po to satisfy those constraints. After reshaping, Po becomes Lo, and the 
poles /zeros /gain added to compensate Po forms the controller, G 0 as can be 
depicted from this relationship: 

Lo = Po * G c 

To satisfy those constraints, each frequency coj on L 0 should be kept on and above it 

corresponding B(ja)i); L 0 should also not penetrate the U contour. Figure B-5 shows 
the U contour, performance bounds and the optimal Lo of an example problem. 

5) After reshaping(compensation), the system is guaranteed robust over the full range of 
plant uncertainty, i.e. 6-r(jwj) <; SR(jwj). However, the system may not have met the 
performance specifications completely yet(see Figure B-6.) A prefilter may be 
required to reposition the compensated system to fully meet the specification. 

















